THREE-DIMENSIONAL UNSTEADY SIMULATION OF AERODYNAMICS AND 
HEAT TRANSFER IN A MODERN HIGH PRESSURE TURBINE STAGE 


Unsteady 3-D RANS simulations have been performed on a highly loaded transonic turbine stage and 
results are compared to steady calculations as well as to experiment. A low Reynolds number k-£ 
turbulence model is employed to provide closure for the RANS system. A phase-lag boundary condition 
is used in the tangential direction. This allows the unsteady simulation to be performed by using only one 
blade from each of the two rows. The objective of this work is to study the effect of unsteadiness on rotor 
heat transfer and to glean any insight into unsteady flow physics. The role of the stator wake passing on 
the pressure distribution at the leading edge is also studied. The simulated heat transfer and pressure 
results agreed favorably with experiment. The time-averaged heat transfer predicted by the unsteady 
simulation is higher than the heat transfer predicted by the steady simulation everywhere except at the 
leading edge. The shock structure formed due to stator-rotor interaction was analyzed. Heat transfer and 
pressure at the hub and casing were also studied. Thermal segregation was observed that leads to the 
heat transfer patterns predicted by steady and unsteady simulations to be different. 
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NASA Subsonic Transport System Level Metrics 

.... technology for dramatically improving noise , emissions, & performance 


CORNERS OF THE 
TRADE SPACE 

N+1 (2015)*** 
Technology Benefits 
Relative to a 
Single Aisle Reference 
Configuration 

N+2 (2020)*** 
Technology Benefits 
Relative to a 

Large Twin Aisle Reference 
Configuration 

N+3 (2025)*** 
Technology Benefits 

Noise 

(cum below Stage 4) 

- 32 dB 

- 42 dB 

-71 dB 

LTO NOx Emissions 
(below CAEP 6) 

-60% 

-75% 

better than -75% 

Performance: 
Aircraft Fuel Burn 

-33%** 

-40%** 

better than -70% 

Performance: 
Field Length 

-33% 

-50% 

exploit metroplex* concepts 


*** Technology Readiness Level for key technologies = 4-6 

** Additional gains may be possible through operational improvements 

* Concepts that enable optimal use of runways at multiple airports within the metropolitan areas 


SFW Approach 

- Conduct Discipline-based Foundational Research 

- Investigate Advanced Multi-Discipline Based Concepts and Technologies 

- Reduce Uncertainty in Multi-Disciplinary Design and Analysis Tools and Processes 

- Enable Major Changes in Engine Cycle/Airframe Configurations 
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About this work... 


• 3-D URANS simulations performed on highly loaded transonic turbine stage 
using TURBO (Chen et al.) 

• Results are compared to steady calculations as well as experiment (Tallman, 
Haldemann et al. at OSU GTL). 

• Effect of unsteadiness studied 

- shock structure 

- rotor heat flux 

- hub and casing heat flux 

- thermal segregation 
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Background 


• HPT flow is unsteady due to wake passage and shock-wake interactions 

• Stagnation point on rotor moves away from LE 

• Shock moves from rotor crown to leading edge as wake passes (Denos et al., 
Paniagua et al.) 

• Thermal segregation could occur (Shang and Epstein, Ameri et al., Kerrebrock and 
Mikolajczak) 
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Segregation 




Lower temperature in wake directs cooler gas to suction side 
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TURBO 


• Upwind Roe scheme with Newton sub iterations 
-No artificial dissipation 

• Fully parallelized to use MPI 

• Only need one blade per row using Phase lag 

• Phase lag - ideal for single-stage simulation (Van Zante et al.) 

• Uses low Re k-s turbulence model 

• Heat transfer simulation made possible by incorporating isothermal BC 


Phase lag 

Uses blade count of 
neighboring blade row to 
determine frequency 


Fundamental Aeronautics Program 
Subsonic Fixed Wing Project 


6 



The Grid 



Rotor grid. 


Rotor Features 

- 2.1% tip clearance 

- blade speed - 9000 rpm 
-Re~3x10 6 /m 


Convergence 
determined using mass 
flow and surface heat 
transfer 


- 38 Stators, 72 rotor blades 

- ~2.5M cells (very fine grid) 

- y+<1 (at first point off wall) 

- Coarser grids have shown satisfactory results 
(e.g. Green et al.) 



Pressure 

Exit 

(profile 

specified) 


Isentropi 

inlet 


Time 

shift 


Sliding 

interfaces 


Relative stator-rotor positioning for unsteady case and 
boundary conditions 
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Simulations 


Isentropic 

inlet 



Sliding 

interfaces 


Pressure Exit 

(profile 

specified) 


Time shift 


Steady 

- Used circumferentially averaged vane exit 
total pressure and temperature profiles as 
rotor inlet profiles. 

- Periodic BC used in tangential direction 
Unsteady 

- Stator inlet total pressure and temperature, 
rotor exit static pressure specified 

- Phase lag BC in tangential direction 

- 50 steps per period stored 
Sliding interface BC at stator-rotor interface 
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RESULTS 
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Shock Function at 15% span 


SF = V ■ Vp 


•Boundaries of 
red regions are 
shocks 

(Large pressure 
gradient is in 
direction of flow 
velocity) 



| 1.265 

®- 1.267 
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Shock Function at Mid-Span 
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Vane trailing edge 
shock sweeps 
across crown and LE 
of rotor 

^ Minimum 

unsteadiness on 
suction side of rotor 

1 .267 
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Comparison with Schlieren 






De la Loma, Paniagua, Verrastro, Adami, 
GT2007-27 101 



Dotted lines are 
reflected shocks 
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Shock Function at 90% Span 


Vane trailing edge 
shock interacts with 
rotor trailing edge 
shock 
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Mid-span Pressure and Mach Number 




Pressure 


Mach no. 
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Pressure 



A 


a) 15% 


P/P 


ref 


0.7 
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Steady CFD 

Time-Averaged CFD 

-♦-Unsteady Envelope Max 

Unsteady Envelope Min 

• Experiment 
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s=o 


S (non dimensional distance along blade profile) 



S= -1 S=1 


Unsteady envelope widest near LE, 
narrows towards TE 


Rotor 
trailing 
edge shock 
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Pressure Profiles 


b) 50% 


P/Pref 


c) 90% 



- 1.00 - 0.80 - 0.60 - 0.40 - 0.20 0.00 0.20 0.40 0.60 0.80 1.00 

S (non dimensional distance along blade profile) 


0.70 


0.60 
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0.40 


0.30 


0.20 


0.10 


Steady CFD 

Time-Averaged CFD 

—•—Unsteady Envelope Max 

Unsteady Envelope Min 

Experiment 


0.6 


O.S 


0.4 


0.2 


0.1 


0.00 

- 1.00 - 0.80 - 0.60 - 0.40 0.20 0.00 0.20 0.40 0.60 0.80 1.00 


0.3 
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Stag. Point 
moved to 
pressure 
side 



- No shock near 
tip 

- Thin envelope 
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Stanton Number 


A 

a) 15% 

St 



0.050 
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0,040 
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Steady CFD 

Time -Averaged CF O 

unsteady Envelope Man 
-• Unite a dy E nve lope M-in 
* Experiment 
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S (non dimensional distance along blade profile) 
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Stanton Number 


b) 50% 

St 

c) 90% 
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Streamlines 



Streamlines of relative velocity over suction side of rotor blade with rotor blade 
showing Stanton number contours. 
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Surface Heat Flux 



a) Steady b) Time-averaged 

Comparison between steady and time-averaged Stanton number distribution on rotor 
blade pressure side. 
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Snapshots of Unsteady Heat Flux 
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In the Tip Gap 



Stanton number on blade 


0 = tan l (— ) 

w 

Isosurface :x — z tan(0) = c 

V plane = « Sill(<?) + WCOS(<9) 


Vpiane on Iso planes of c 


Vpiane is positive (yellow through red 
colors) if flow is separated. This is only 
valid in the tip gap, where flow 
separation can be measured in the z-x 
plane 
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In the Tip Gap - Plane 1 


Expands to pressure lower than suction side 
and then goes through series of 
compressions and expansions 
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In the Tip Gap - Plane 1 - Unsteady 



Shock function 


tip_shock3.avi 
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Casing Heat Transfer 


1 


0.040 


0,000 


Steady 


Time- 

averaged 




Corresponds to 
separation 
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Hub Heat Transfer 


St —St 

steady u L time-averaged 


St 


time-averaged 


I 100.000 

33.333 

- 33.333 

- 100.000 



Segregation 

effect 


Blue: Time- 
averaged 
heat 
transfer 
higher 


Percent difference between steady and time-averaged Stanton number on rotor hub. 
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Conclusions 


• Over most of blade surface, steady simulation is accurate 

• Thermal wake causes unsteady heat transfer over most of the blade to 
be higher than steady heat transfer except at leading edge. 

• At the leading edge the effect of unsteadiness is most prominent 

• Thermal redistribution was observed at the hub and on the blade 
surface 

• Pressure and heat transfer distribution over blade is highly 3D 
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Backup slides 
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Time averaged P and Shock @ 50% span 




— 0.873 
0.222 
*- 0.429 


-1 .079 
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The phase lag boundary condition for more than two blade row 



• In this example, adding the IGV wakes creates circumferential non-uniformities 
at the entrance to stator 1 . 

• At ‘B’ the phase lag boundary condition will apply the time history of ‘A’ with a 
phase shift but not the necessary change in the mean. 

• This results in a spatial filtering of information for stator-stator (and rotor-rotor) 
interactions. 
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Pressure and Temperature Profile at Interface 
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Stanton No. Derivation 


St = 


h 

Cp-p-V 


since h = T . and knowing lii = p - A - V 

‘wall- ‘ ref ° r 


St = 


tj wall 


(f ) ■ c P ■ [T wall - T ref ] 


Pr — 


Inwall 't- P 


St = 


1 


P'V Cp T W all - T re f 


-k ■ Sr \waii knowing: 


k 

n - Pr k 
Mu 'all ~ Cp 

t^wall Pwall ' V\vaU 


multiply and divide by j u wa ii and simplify to obtain: 

dT 


dn 


wall Pwall I 


Twall Tre f p ■ V Pr 
An expression in nondimensional terms: 


Sr = - 


dT 

dn 

wall 


Inwall 


‘ 1 ‘ 


Pref 

Twall 

— T rc f 


p • V 


Pr 


_ Preface 0 _ 


Suosomc riAeu winy rrujeci 


Rotor Hub Surface Heat Transfer - Steady 


TURBO predictions Stanton No. Comparisons 



Fundamental Aeronautics Program 
Subsonic Fixed Wing Project 


Rotor Tip Heat Transfer 


TURBO predictions 



0.012 

0.010 

0.008 

0.006 

0.004 

0.002 

0.000 
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Steady 


Stanton No. Comparisons 


Stanton Number vs. Probe Number 
Rotor Tip 
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♦ Experiment #5 
■ Experiment #6 
ATallman- CFD 

• Turbo - CFD 

1 — 






0 1 2 3 4 5 

Probe Number 


Rotor Casing Surface Heat Transfer - Steady 


No data available for 

TURBO predictions comparison 



3 2 1 


■ 0.04 

Stanton Mo. 



Stanton Number vs. Probe Number 



Probe Number 
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Statistics - Steady 


Vane Surface Pressure 
Vane Surface Heat Transfer 
Blade Surface Pressure 
Blade Surface Heat Transfer 

Rotor Hub Surface Heat 
Transfer 


Good 

Good 

Good 

Fair 

13.4% difference 


Good 

Good 

Good 

16.3% difference (max: 
25%) 

13.1% greater 


Rotor Tip Heat Transfer 


10.5% difference 15.5% greater 


Stator 30,000 60 Hours 

Rotor 40,000 120 Hours 
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